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Abstract —We address the max-min fairness design problem 
for a MU-MISO system with partial Channel State Information 
(CSI) at the Base Station (BS), consisting of an imperfect channel 
estimate and statistical knowledge of the estimation error, and 
perfect CSI at the receivers. The objective is to maximize the 
minimum Average Rate (AR) among users subject to a transmit 
power constraint. An unconventional transmission scheme is 
adopted where the Base Station (BS) transmits a common 
message in addition to the conventional private messages. In 
situations where the CSIT is not accurate enough to perform 
interference nulling, individual rates are assisted by allocating 
parts of the common message to different users according to 
their needs. The AR problem is transformed into an augmented 
Average Weighted Mean Square Error (AWMSE) problem, solved 
using Alternating Optimization (AO). The benefits of incorporat¬ 
ing the common message are demonstrated through simulations. 

Index Terms —MISO-JMB, max-min fairness, AWMSE. 

I. Introduction 

We consider a Multi-User (MU) Multiple Input Single 
Output (MISO) system where a Base Station (BS) equipped 
with N t antennas serves K < N t single-antenna receivers. For 
a given transmit power P t , the objective is to design precoders 
(or beamformers) that maximize the minimum rate among 
users; a criteria known in literature as max-min fairness JT)- 
0. In practical scenarios, perfect Channel State Information 
at the Transmitter (CSIT) is unguaranteed due to the varying 
nature of wireless channels and the limited resources available 
for CSIT acquisition. This prompted a number of robust 
fairness-based designs 0-®. In this paper, we consider the 
case where the BS obtains an imperfect estimate of the channel 
(through feedback or training) and statistical properties of 
the estimation error, while users estimate and track their 
channels accurately. The problem is formulated in terms of 
the Average Rate (AR) which captures the overall performance 
w.r.t channel estimation errors. Contrary to our previous work 
in (6l which considers a similar scenario, a rather unorthodox 
transmission scheme is considered in this paper. In particular, 
the BS transmits a common message in addition to the 
conventionally transmitted private messages. This scheme is 
termed Joint Multicasting and Broadcasting (JMB). 

JMB was employed in 0 to boost the achievable Degrees 
of Freedom (DoF) under estimation errors that decay with 
increased Signal to Noise Ratio (SNR). Assuming fixed noise 
variances at the receivers, the decaying power of the estimation 


error scales as 0(P t _a ), where a > 0 is an exponent that 
represents the CSIT quality. This power relationship can be 
used to express a variety of imperfect CSIT scenarios f8j]. For 
example, a = 0 represents the case where the error power 
is fixed over the whole SNR regime, e.g. constant number 
of feedback bits. On the other hand, a = oo corresponds 
to perfect CSIT. In DoF analysis, it is customary to truncate 
the exponent such that a < 1, where a = 1 is considered 
perfect from a DoF perspective as the error impact becomes 
negligible for P t —> oo. Under this assumption, conventional 
MU transmission that employs Zero Forcing (ZF) Beamform¬ 
ing (BF) achieves a sum DoF of Ka. Interestingly, the same 
sum DoF can be achieved using only part of the power given 
as P['. This can be explained as follows: in the presence 
of CSIT uncertainty, perfect interference nulling cannot be 
achieved, and any DoF gain achieved by increasing the power 
above a certain threshold (i.e. P“) is nullified by the increased 
interference. Nevertheless, the remaining part (i.e. P t — ) 

can be used to transmit a common message which achieves a 
DoF of 1—a by decoding it while treating the private messages 
as noise. As the common message is decodable by all users, 
it can be cancelled from their received signals, and the DoF 
achieved by the private messages is unaffected. 

Contribution and Organization: We employ JMB to max¬ 
imize the minimum AR among users, where a given user’s 
AR consists of a private AR, in addition to a portion of 
the common AR. This is achieved by jointly optimizing the 
precoders and the partition factors which determine how much 
of the common message is allocated to each user. This non- 
convex problem is transformed into an equivalent augmented 
Average Weighted Mean Square Error (AWMSE) problem, 
which is solved using Alternating Optimization (AO). Simula¬ 
tion results show the superiority of the proposed scheme over 
conventional MU transmission, which does not incorporate the 
common message. For example, while conventional transmis¬ 
sion with a = 0 yields flat rates in the high SNR regime (6J, 
JMB employs the multicast part to achieve a non-vanishing 
DoF. To the best of our knowledge, multicast assisted robust 
fairness-based designs have not appeared in literature. 

The rest of the paper is organized as follows: the system 
model and problem formulation are given in Section UD The 
equivalent AWMSE problem is formulated in Section |III1 The 
AO algorithm is proposed in Section [IV] Simulation results 
are given in Section [VI and Section [VI] concludes the paper. 


Notation : Boldface uppercase letters denote matrices, bold¬ 
face lowercase letters denote column vectors and standard 
letters denote scalars. (-) T and (-) H denote transpose and 
conjugate-transpose (Hermitian) operators, respectively, tr(-) 
and || • || are the trace and Euclidian norm operators, respec¬ 
tively. E.„ {■} is the expectation w.r.t the random variable x. 

II. System Model and Problem Formulation 

For the MISO setup described in the previous section, 
the vector of complex data symbols is given as s = 
[s c , Sl ,...,s K } T £ C K+1 , where s c is the common symbol 
decodable by all users, Si is a private symbol intended solely 
for the ith user, * £ 1C and 1C = {1,..., K }. Entries of s have 
zero-means, unity powers and are mutually uncorrelated such 
that E{ss h } = I. s is linearly precoded as 

K 

X = PcS c + ^p lSl Cl) 

2—1 

where x £ C Nt is the transmit vector. p c £ C Nt and 
p, £ C Nt are precoders for the common symbol and the ith 
private symbol respectively, from which P = [p c , Pi, ■ ■ ■, Pa] 
is composed. The BS transmit power constraint is given as 
E{x H x}=tr(PP fl ) < 1\. The A'th received signal is given as 
y k = hf x + n k (2) 

where h/ ; . £ C Nt is the channel impulse response vector 
between the BS and the A'th user, from which the composite 
channel is defined as H = [hi,..., Iir-]. n k ~ CA/”(0, er 2 fc ) 
is the AWGN. Throughout the paper, it is assumed that noise 
variances are equal across all users, i.e. er 2 fc = cr 2 , Vk £ 1C. 

A. CSIT Uncertainty 

Each of the K links exhibits independent fading, and 
remains almost constant over a frame of symbols, enabling 
users to estimate their channel vectors with high accuracy. 
On the other hand, CSIT experiences uncertainty arising from 
limited feedback, delays or mismatches. H is written as a sum 
of the transmitter-side channel estimate H = [hi,..., Iir-] 
and the channel estimation error H = [h |..... li/ N -], such 
that H = H + H. The CSIT consists of H, in addition to 
some statistical knowledge of H. Particularly, the BS knows 
the probability distribution of the actual channel given the 
available estimate, i.e. / H |h(H) = /g(H — H). 

B. MSE and MMSE 

The fcth user obtains an estimate of the common symbol by 
applying a scalar equalizer g c>k (h k ) to ([2} such that s Ci fc = 
<?c,fc(hfc)t/fc. Assuming that the common symbol is successfully 
decoded by all users, the common symbol’s receive signal 
part is reconstructed and cancelled from y k . This improves 
the detectability of s k , which is then estimated by applying 
g k ( hfc) such that s k = g k (h k )(y k - hf p c s Cjfe ). The notations 
3 c,fc(hfc) and g k (h k ) are used to emphasise the dependencies 
on the actual channel, as each user is assumed to have perfect 
knowledge of its own channel vector, (h*,) is omitted for 
brevity unless special emphasis is necessary. This is also used 


with other variables that depend on the actual channel. For the 
fcth user, the MSEs defined as s Ctk = E Sjnfc {|s" C) fc — s c | 2 } and 
£k = E S! „ fc {|s fe - s fc | 2 } are given as 

£c,fc(h fe ) = \g c ,k\ 2 T Cyk - 25R{g c ^hf p c } + 1 (3a) 

e k (h k ) = \g k \ 2 T k - 23?{ fffe hf p fc } + 1 (3b) 

where T C)fe = |pf h fc | 2 + T fe and T k = |pfh fc | 2 + cr 2 . 

Optimum g c , k and g k are obtained from gf’} =0 and =0, 
yielding the well-known Minimum MSE (MMSE) equalizers: 


ffc M , MSE (h fe )=pfh fc T- fc 1 and ff f MSE (h fc )=pfh fc T- 1 . (4) 
Substituting (|4} into (O, the fcth user’s MMSEs are given as 

£c M f SE (h k ) = TjE Ctk and ef MSE (h k )=T^E k (5) 

where = T Cjfe - |pf h fc | 2 = T k and E k =T k - |pf h k \ 2 . 


C. Achievable Rate and Average Rate 

The MMSE and the Signal to Interference plus Noise Ratio 
(SINR) are related such that 7 Cjfc = (1 — e^ 1 f ISE )/e^ I f ISE and 
= (1 — ef MSE )/ef MSE , where y Cyk and y k are the fcth 
user’s SINRs. Therefore, the fcth user’s maximum achievable 
common rate and private rate are written as R c k (h k ) = 
— log 2 (eff ISE ) and R k (h k ) = — log 2 (£f MSE ), respectively. 
The common message is transmitted at a common rate defined 
as R c = min j{Rc,j}jLi, which ensures that it is decodable by 
all users. As the multicast part is used to boost the individual 
rates achieved by users, the fcth user is allocated a fraction c k 
of the common message, where c k > 0 and °k = 1- 

The corresponding portion of the common rate is given as 
c k R c . In this case, the fcth user’s total achievable rate is given 
as R k + c k R c . In a block-based transmission, this can be 
achieved by dedicating different parts of the common data 
block to different users, where any fraction c k can be achieved 
for sufficiently long blocks. To achieve fairness among users, 
the objective would be to design P and c = [ci, ..., cat] t 
that maximize the minimum total rate among users defined as 
R = minfc {R k + c k R c } k _ . However, rates depend on the 
actual channel, and hence cannot be considered as optimization 
metrics at the BS. Alternatively, we consider the Average Rates 
(ARs) defined as: E h |} and E^^ {R k }, where 
averaging is taken over the CSIT error. In the following, 
Eh fc |h fc {’} simply referred to as E{-}. 


D. Sample Average Function 

In order to formulate a deterministic problem, the stochas¬ 
tic ARs are approximate by corresponding Sample Average 
Functions (SAFs) obtained by taking the ensemble aver¬ 
age over a sample of M independent identically distributed 
(i.i.d) realizations drawn from the distribution / H |g- The 
sample is defined as H_vf = |H ( m ) | m 

£ A4 }, where 


is the mth Monte-Carlo real- 


H M a [ h (" 

ization, and M = {1..... A/}. The SAFs are given as: 
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TDK 

m—1 ,k 


Am) 


and R k 

? (m) 
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where 


and R { j^ n) = R k ( h} m ^) are the rates 


associated with the realization h}" *. In the following, the 



superscript (to) is used to indicate the association of variables 
with the ruth Monte-Carlo realization. It should be noted that 
P is fixed over the M realizations of the rates, which follows 
from the definition of the ARs. This also reflects the fact that 
P is optimized at the BS using partial CS1 knowledge. 

Assumption 1. In the following, we assume that || hfc | p Pt > 0 
with probability 1, Vi £ K. 

Alternatively, we can say that SNR = P t /cr(; can only grow 
finitely large, and channel gains are finite. Assumption|T|yields 
£ ™ SE i £ f MSE > 0 with probability 1, as the presence of a 
nonzero noise variance dictates that E ck ,E k > 0. This also 
implies that rates are finite, and by the strong law of large 
numbers we can write 

R c ,k— A ]im R { c M k ] = E{P C , fc }, almost surely (6a) 

R k = lim R = E{Pfc}, almost surely (6b) 

M—too 

where R ck and Rk are the approximated ARs for a sufficiently 
large M, which will be simply referred to as the ARs. More¬ 
over, the common AR is defined as R c = niin ; { R (J }jLj ■ 

E. Problem Formulation 

The minimum total AR among the K users is defined as 
R = min k { A/,-+c/ r /?c} ; . r The objective is to design P and 
c such that R is maximized. This problem is formulated as 


1Z : max R 

(7a) 

R,R c ,T?,c 


s.t. R k + c k R c > R, Mk £ /C 

(7b) 

Rc,k ^ Rc, vfc g /c 

(7c) 

Cfc >0, Vfc G /c 

(7d) 

K 


H cfc = 1 

(7e) 

k= 1 


tr(PP H ) < P t 

(7f) 


where the constraints in (TTbt and (ITcl) are introduced to elim¬ 
inate the potential non-smoothness arising from the pointwise 
minimizations in R c and R. 7Z is non-convex and appears to 
be very challenging to solve. 

III. Equivalent AWMSE Optimization Problem 

In this section, 1Z is transformed into a more tractable 
equivalent problem by exploiting the relationship between 
rates and augmented WMSEs, an approach inspired from 0, 
0. We start by introducing the main components used to 
construct the equivalent problem, i.e. the augmented WMSEs: 

Cc,fc Pc,fc) — t£ C ,A' (Ilfc)^c,fc(hfc) log 2 (hfc)) (8a) 

£k(h k ,g k ,u k ) = u k (h k )£ k (h k ) log 2 (ufc(h fc )) (8b) 

where u C! k{iik) > 0 and u k (h k ) > 0 are weights associated 
with the /, th user’s MSEs. The dependencies of E, c .k and £ k on 
different variables are highlighted in (0 for their significance 
in the following analysis, where we establish the following 
WMSE-Rate relationship: 


min = 1 - P c , fc and min £ k = 1 - R k . (9) 


This can be shown as follows. From d .^ c ' k = 0 and = 0, 

og c ,k ogk 

the optimum equalizers are obtained as g* k = p^ ISE and 
9k = Substituting this back into 0, we obtain the 

augmented WMMSEs written as 


&“fc MSE ( h fc, u c,k) = u Ctk e^ SE - log 2 (u Cjfc ) (10a) 
£f MSE (h k ,u k ) = u k ef MSE - log 2 (u fc ). (10b) 


gcMMSE ^>MMSE 

Furthermore, from — = 0 and -%- = 0, we obtain 

ou c ,fc ouk 

the optimum MMSE wights: u* k = u^ I ^ ISE = (e^ [SE ) 1 
and u k = u™ MSE = (efc IMSE ) '. Substituting this back into 
([Tol l yields the relationship in 0. It is evident from 0 that 
the MMSE weights are dependent on the channel. 

The equivalent problem is formulated using the augmented 
AWMSEs defined as: E{£ Cj fc} and E{£ k }. Before we proceed, 
the augmented AWMSEs are approximated as: 

M i M 

& = jjJ2 & and li M) = mJ2 £l ro) > where 


m= 1 


& - >SS “S’) -a «r=& (it 


m— 1 

(m) A. c ( u {m) 
k i 


9 k 


(m) (m)\ 


correspond to the ruth realization of the augmented WM¬ 
SEs, which depend on the ruth realization of the equalizers: 
9^k = p c ,fc(h fe m) ) and g[ m) = g k { h^ m) ), and the weights: 
u^ k = u c ,fc(hfc m) ) and = zifc(hj. m) ). For compactness, 
we define the set of equalizers associated with the M real¬ 
izations and the K users as: G = {g c ,fc,gfc | k G /C}, where 
gc,fc = I rn G M] and g fc = {g ( k m) \ m G M). In a 

similar manner, we define: U = {u c fc,Ufc | k G /C}, where 
u c ,fc = {u^™k | to G M} and u fc = {u^ \ m G X}. The 
approximated augmented AWMSEs for a sufficiently large M 
are defined as ^ Cjk = lim. M -koo^P and £ k = limM->oo^ M) , 
which will be simply referred to as the AWMSEs. The same 
approach used to prove 0 can be employed to show that 

min Cc.fe = 1 - Rc,k and min £ k = 1 - R k . (11) 

U c ,fc,gc,fe ?Sfe 


where optimality conditions are checked separately for each 
realization. The sets of optimum MMSE equalizers associated 
with (fTTb are defined as g^ 13 ® = | m g yVf} 

and g^ MSE = |g^ MSE ( m ) | m g _A/fj. j n j-jjg same man¬ 
ner, the sets of optimum MMSE weights are defined as 

u MMSE (m) A { u MMSE (ra) | m g an( j u MMSE( m ) A 

{uf MSE(m) 1 171 ^ • F° r t ^ le ^ users, the MMSE solution 

is composed as G MMSE = {g^ 1 ^ ISE , g^ MSE \ k G and 
IJMMSE a { u mmse , u MMSE | k € k ). 


A. Augmented AWMSE Minimization 

Prompted by the relationship in ifTTb . the augmented 
AWMSE minimization problem is formulated as 






min £ 

(12a) 

£,£ c ,P,c,U,G 


Cfc + c — l) < £, Vfc G /C 

(12b) 

fc ,k < Ccj Vfc G /C 

(12c) 

Cfc > 0, Vk G K, 

(12d) 

K 


Y, Ck = 1 

(12e) 

k =1 


tr(PP ff ) < P t 

(12f) 


where (£, £ c ) are auxiliary variables. The equivalence between 
A under the MMSE solution, and 7Z is demonstrated as fol¬ 
lows. Substituting (g^ MSE , u]4 MSE ) and (g^™ SE , u^™ SE ) 
into ( 1 1 2bb and ( 1 1 2cb respectively, and rearranging, we obtain 
Rk+c k (l-£ c ) >1-1 V/cG/C (13a) 

R c ,k >1 - VA G 1C. (13b) 

It can be seen that (1 1 3ab and ( 1 13 IT are equivalent to (TTbb and 
(I7cl respectively, where 

R= l-Z and R c = 1-&. (14) 

Moreover, the monotonicity of the relationships in (fl4l > im¬ 
plies that the two problems are equivalent under the MMSE 
solution, i.e. maximizing R is equivalent to minimizing £. This 
relationship is exploited in the following section. 

IV. Alternating Optimization Algorithm 

Although A is non-convex in the joint set of optimization 
variables, it is convex in each of the blocks G, U, c and 
P. We propose an AO algorithm that exploits this block-wise 
convexity. Each iteration of the algorithm consists of three 
steps: 1) updating G and U, 2) updating c, 3) updating P. 

A. Updating the Equalizers and Weights 

In nth iteration of the AO algorithm, the equalizers 
and weights are updated as: G = QMMSE^p) anc j 
U = U MMSE (P), where P is the precoding matrix ob¬ 
tained from (n — l)th iteration. For the resulting point 
(P,G mmse (P),U mmse (P)), the updated AWMSEs in 
are denoted by f.j. and f k . Moreover, the common AWMSE 
is given as £ c = lniiij {£ CJ }f —\. Finally, the cost function in 
( fl2T> is reduced by minimizing the individual AWMSEs. 

B. Updating Partition Coefficients 

Using the updated AWMSEs from the previous subsection, 
and employing the relationships in ( ITTI i and ( 1 1 4k the problem 
of optimizing c is formulated in terms of the ARs as 


:min — R 

(15a) 

R, c 


Rk + CkRc ^ R-, V/c G /C 

(15b) 

K 


Z )Ck = 1 

(15c) 

k =1 


Ck > 0, VA: G /C 

(15d) 


where R k = 1 — 0,: and R c = 1 — £ c . While optimizing 
c does not influence the private ARs or the common AR, 


it redistributes the common AR among users in a way that 
further reduces the cost function of A. A solution for problem 
Ac is developed by examining the KKT optimality conditions. 
The Lagrangian of (IT5l > is given as 

k .. .. K K 

L — R -(- y ' pki^R R k Cfci? c ) T A ( y ' J Ck l) y \ r'kttk 

k =1 k-1 k =1 

where pt = {pk \ k G /C}, A, and v = {vk \ k G K.} are the 
multipliers associated with the constraints in (U5ab . (1 1 5bb . and 
(fl5ct . respectively. The KKT conditions can be written as 

d p ^ 

=0 y; Pk = 1 (16a) 

k= 1 

77 — = 0 =► ~pkR c + X - v k = 0, Vfc G K (16b) 
oc k 

Pk > 0, Pk(R - Rk- c k R c ) =0, Vfc G K, (16c) 

v k > 0, v k Ck = 0, V/c G 1C. (16d) 

For R c = 0, optimizing c is irrelevant. On the other hand, 

R c > 0 implies that A > 0. This is concluded from (116bl) 

which can be rewritten as A = Vk+p k Rc , (U 6 al) which dictates 
that at least one element in // is nonzero, and the positivity of 
Vk front (1 1 6db . Now we consider c k under three cases: 

1) For the first case, we assume R < R k . In this case, (116cl) 
dictates that p k = 0, as R — R k — CkR c = 0 is impossible due 
to the positivity of c k R c - Moreover, we have A = v k > 0 from 
(fl 6 bt . This requires c k = 0 as seen from (1 1 6db . In other words, 
if a given user’s private AR is higher than the optimum R, it 
will not be allocated a fraction of the common AR. Moreover, 
the corresponding constraints in (1 1 5bb are inactive. 

2) Next, we consider R = Rk ■ Assuming that c k > 0, then we 
must have pk = 0 to satisfy (II 6 cb . This implies that \ = v k > 0, 
which contradicts (1 1 6db . Therefore, we must have c k = 0 and 
R=Rk, i.e. the corresponding constraints in (1 1 5bb are active. 

3) Finally, we consider R > R k . Assuming that //./,. = 0, 
then we have A = v k > 0 and c k = 0. However, this is 
an infeasible solution as it contradicts (fl5bt . Therefore, we 
must have pk > 0 and R = R k + c k R c - For this case, the 
corresponding constraints in (1 1 5bb are also active. 

From the previous analysis, we can write 

c k =ma,x{0,Rf 1 (R-Rk)}, VkeK, (17) 

which suggests that the optimum c can be calculated using 
a form of water-filling, where R is the water level. Without 
loss of generality, we assume that the private ARs are ordered 
in an ascending manner, i.e. Ri > Rj , V'i, j G K and 
i > j. Assuming that all constraints in (U5bb are active, and 
using the constraint in (fl5ct . the water level is calculated 
as R = K~ 1 (R C + Rk) , from which we obtain 

Ck = Rf 1 ^—Rk), Vfc G 1C. If the /\th user (with the greatest 
private AR) has ck < 0, then the first case of the three applies. 
The A’th user is discarded and R and c are recalculated for 
the remaining K — 1 users. Otherwise, the optimum solution 
is obtained. This is summarized in Algorithm [j] 
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Algorithm 1 Optimizing the partition coefficients, 
l: Initialize: K c <— K + 1 

2: repeat 

3: K c K c — 1 and IC c <— {1,..., A c } 

4: R^K-\R,+^ l R k ) 

5: Ck R c 1 (A — Rk) • Vfc S /C c 

6: C k 0. Vfc S K. \ /C c 

7: until C/f c > 0 


C. Updating the Precoders 

To formulate the problem of updating the precoders, the 
AWMSEs are written in terms of the updated blocks G and U, 
and the block P which is yet to be updated. For this purpose, 
we introduce the AWMMSE-components: T , c .fc, 'I 'k, t c ,k, tk, 
fc.fc, f k, Uc,k, Uk, v Ci k and Vk, which are obtained using the 
updated G and U. In particular, u c j i: and Uk are calculated 
by taking the ensemble averages over the M realizations of 
u™ and u^ n> . The rest are calculated in a similar manner by 
averaging over their corresponding realizations: 


Am) 

,k 

(m) (m) ^ 

= U c,k 9c, k 

and 

Am) 

l k 

(m) (m) ^ 

= u k 9l 

*( m ) 
^ c,k 

~ L c,k lL k ll k 

and 


_ ,(m), (m), (m) 
— L k ll k ll k 

Am) 

L c,k 

— 

— a c,k ll k 9c,k 

and 

f (m) 

l k 

~ a k ll k 9k 

(m) 

V o,k 

= l0 S2 ( U cl } ) 

and 

( m ) 
V k 

= log, • 


Using the AWMMSE-components, the AWMSEs are given as 

K 

Cc ,k = pf’S'c.fePc + pf ^c.fcP i + ,k - 2Sft{f c 5.p c } 

2=1 

H - ,k ,k (18&) 

K 

Cfc = ^pf^kPi + ttltk - 25R{ff p fc } +u k -v k (18b) 

i=l 

from which the problem of optimizing P is formulated as 
,4p : min £ (19a) 

K 

s.t. Y i p?<S’kP i +a 2 n tk-2M{WPk}+Uk-Vk 

2=1 

+c fe (ec - 1) < £, Vfc G K, (19b) 

K 

Pc^c,kPc+'^2pf^fc,kPi+ctli c ,k-2^i{i^kPc} 

i=1 

+ t£c,fc V c ^k A £ c , Vfc C /C (19c) 

tr(PP ff ) < P t . (19d) 

Problem ( IT9t is convex with quadratic constraints, which can 
be solved using off-the-self optimization software that employs 
interior-point methods iflOl . 

D. Alternating Optimization Algorithm 

The AO algorithm is constructed by repeating the steps 
described in the previous subsections until convergence. This 
is summarized in Algorithm [2] where e r determines the 



Fig. 1. Convergence of Algorithm [2] using 1 randomly generated H. 


accuracy of the solution and n max is the maximum number of 
iterations. Initializing P is discussed in Section [VI 

Algorithm 2 Alternating Optimization 


Initialize: n «— 0, «— 0, P 

repeat 

n n + 1, P f- P 
G -f- g mmse (P), U <- u MMSE (P) 

c 4— arg Ac 

Update {(Pc,/;;, (Pfc j Pc,fc i Pfc i f c,k : f k i ; FI k , V c k > Ft A;} 

(P,c (n) ) <- arg Ap and R ("> <- 1 - ^ 
until M n ^ — R( n ~ 1 ' > < eg orn = n max 


fc=i 


Proposition 1. Algorithm [2] monotonically increases the AR 
objective function of problem 7 Z until convergence. 

Proof: From the steps in Algorithm [3 we observe 
the sequence: 1 — < R < R < 1 — £("), where 

R = vam.k{ckR c + Rk} k -i * s t ^ le AR objective function at 
the output of step|4l and R = mink{ckRc+Rk} I ^_ 1 is the AR 
objective function after updating c in step 0 The relationships 
in CD, O and (IT4l) hold, particularly at the output of steps 
|4] and 0 This implies that each iteration increases the AR 
objective function. Moreover, the fact that R is bounded above 
for given P t ensures convergence. However, global optimality 
is not guaranteed due to the non-convexity of 1Z. ■ 

V. Numerical Results 

We consider a system with N t = K = 2 and uncorrelated 
channel fading, where the entries of H have a complex 
Gaussian distribution CM (0,1). The noise variance is fixed as 
cr„ = l, yielding a long-term SNR of Pt. Moreover, Gaussian 
CSIT error is assumed where the entries of H are generated 
according to the distribution CM (0, af). For each realization 
H, a channel estimation error H is drawn from CM (0,<Tg), 
from which the channel estimate is calculated as H = H—H. A 
channel realization H should not be confused with a Monte- 
Carlo realization H' m) . The former is the actual channel which 
is unknown to the BS, and assumed to remain constant for a 
given transmission. On the other hand, the latter is part of 
a sample H_vi employed in the optimization at the BS. The 
size of the sample is set to M = 1000. For a given H, the 




























Fig. 2. MISO-BC and MISO-JMB ERs. K = 2, and a\ = P“ 0 ' 6 . 



Fig. 3. MISO-BC and MISO-JMB ERs. K = 2, and cr^ = 0.063. 

mth Monte-Carlo realization is obtained as H fm) — 
where is drawn from CAf (0, <7%). 

First, we examine the convergence of Algorithm U using 
two different P initializations. The first initialization (ZF-e) 
is taken as the sum-DoF motivated design in (7). The private 
precoders are initialized as p/. = \JP™ /Kpf F , where p| F is 
a normalized ZF-BF vector constructed using the channel esti¬ 
mate H. The common precoder is given as p c = \[Rt / J /'ei, 
where ei is a vector with 1 as the first entry and zeros 
elsewhere. The second initialization (ZF-SVD) retains the ZF- 
BF part and the power allocation. However, the common 
precoder is obtained as the dominant left singular vector of H. 
FigureQ]shows the AR convergence for <r f 2 = P t 0 , and SNRs 
5, 20 and 35 dB. It is evident that the algorithm eventually 
converges to a limit point regardless of the initialization. 
Moreover, better convergence performance is achieved with 
ZF-SVD compared to ZF-e due to the initialization of the 
common precoder. The following results employ ZF-SVD. 

Next, we consider the Ergodic Rate (ER) performance. 
For a given channel estimate H and the corresponding error 
statistics, maximizing the minimum AR yields the solution 
P. Employing P for the channel realization H yields the 
achievable rates defined in 1II-CI from which the achievable 
minimum rate among users (i.e. R) is obtained. The ER is 
defined as En{f?}, which captures the overall achievable 
performance for all possible channel realizations. In the fol¬ 


lowing simulations, the ER is calculated by averaging over 
200 randomly generated channel realizations. Conventional 
fairness-based BC transmission is considered as a baseline, 
which is obtained by discarding the common ARs and c in 
1Z. A decaying estimation error power of a 2 = P t ~°’ 6 is used 
in Figure |2] On the other hand. Figure 0uses a fixed estimation 
error power = (10 to) " = 0.063, which is equivalent to 

the CSIT quality obtained at SNR 20 dB in the previous case. 
It is evident that users benefit from the incorporation of the 
common message over the entire SNR range. The significance 
of this utility grows with increased SNR. For example. Figure 
[2] is a manifestation of DoF gains translating into rate gains at 
finitely high SNRs. Moreover, Figure [3] demonstrates JMB’s 
virtue of achieving fairness with non-vanishing growth in a 
scenario where conventional transmission hits a performance 
ceiling. 

VI. Conclusion 

A multicast assisted scheme was proposed to achieve max- 
min fairness in a MU-MISO system with partial CSIT and 
perfect CSIR. Precoders are designed and the common mes¬ 
sage rate is divided among users, such that the minimum AR 
among users is maximized. This problem is transformed into 
an AWMSE problem solved using a converging AO algorithm. 
Simulation results show that the proposed scheme achieves 
significant rate gains over conventional transmission, which 
does not incorporate the multicast part. 
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